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1. Introduction 

Critical phenomena in extended dynamical systems attracted considerable interest 
recently. In particular, coupled map lattices were introduced as paradigmatic models 
of nonequilibrium systems undergoing (in the large-size limit) a second-order phase 
transition between two chaotic states [|l], §]. In those models such a transition is 
associated with a breaking of the Ising symmetry, and is thus expected to share the 
same static critical exponents of the Ising model itself. However, numerical simulations 
0] show that this is not the case. More precisely, the critical exponent of the correlation 
length for the model of Ref. turned out to be significatively different from the Ising 
value. Moreover, the latter are recovered as soon as the updating rule is changed from 
parallel to sequential, thus indicating that some features of the microscopic dynamics 
may be relevant for the critical behaviour. Besides of those facts, it is not completely 
clear why a phase transition behaviour emerges at all and why the coupling is effectively 
"ferromagnetic" . 

For all the above reasons, even the simplest analytical approach, namely the mean- 
field like approximation, is worth to be investigated. Actually, it is not clear to which 
extent it can be applied in nonequilibrium cases and which properties it may share with 
its equilibrium counterpart. Moreover, it is still questionable that such an approach is 
able to reproduce the qualitative behaviour of coupled map lattices in general. 

We will discuss such problems for the coupled map lattice introduced in Ref. 0. 
Its equations of motion read as 

xUi = (1 - 9d)fixl:^) + 9 E ' (1) 

where x^^'^ is the field variable at time n and the index u enumerates the sites on a square 
lattice. Here, Uu denotes the set of the d neighbours to which the z/th site is coupled 
and g is the coupling constant. Rather than considering the original isotropic nearest- 
neighbour coupling on a square lattice we will refer to its modification presented in 
which couples four instead of five sites. This choice has the advantage of simplifying the 
algebraic calculations and offers the possibility to compare our results with the numerical 
simulations reported in |^. Since we employ some mean-field theory, our approach 
applies regardless of the special type of spatial geometry and boundary conditions do 
not play a significant role. 

The local map, which is defined on the interval [—1, 1] 

(-3x-2 for -1 < a; < -1/3 

f{x) = I 3x for -1/3 < a; < 1/3 (2) 

[-3x + 2 forl/3<x<l 

obeys an Ising like symmetry f{x) = —f{—x). In the supercritical region g > Qc the 
symmetry is spontaneously broken and the coupled map lattice exhibits two ordered 
phases, characterised by a non-vanishing value of the "magnetisation" {Y.u^^'^^)- 
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2. Mean— field approach 

Our starting point is the evolution equation for the single-site probability density p{x) 
(cf. e. g. 01), as usually obtained by projection of the Perron-Frobenius equation for the 
full probability density. The single-site density is expressed in terms of the four-sites 
joint probability density p*^^-* as 

pn+iix) = J dyodyidy2dy36 [x - T{yo, yi, 2/2, 2/3)] Pn\yo^ 2/i, 2/2, 2/3) • (3) 
Here we have introduced the abbreviation 

nyo, 2/1, 2/2, 2/3) = (1 - Sg)f ivo) + 9 [/(2/1) + /(2/2) + /(2/3)] . (4) 

The function T is invariant under all permutations of the variables ?/i,?/2,2/3- The 
mean-field Perron-Frobenius equation [^] is as usual obtained by neglecting multiple 
correlations, namely by letting 

Pi^H2/0,2/l,2/2,2/3) Pn(2/o)Pn(2/l)Pn(2/2)Pn(2/3) • (5) 

Then eq.(^ becomes a nonlinear integral equation 

p„+i(x) = / dyodyidy2dy3S [x - T{yo, yi, y2, ys)] Pnivo) pnivi) Pn{y2) Pnivs) ■ (6) 



It should be stressed that the approximation (^) leads to different results if one chooses 
different coordinate systems in the full phase space of the map lattice. In fact, neglecting 
correlations has a different meaning in different coordinate systems. Hence, similar to 
statistical mechanics, it makes no sense to speak about the mean-field approximation. 
The formulation chosen here seems to be quite appropriate for the analytical calculation. 

Since we are interested in stationary properties we look for the fixed point solution 
Pn = P*- In the paramagnetic region, i. e. for sufficiently small coupling, the solution 
is symmetric p*(x) = p*(— x) and does not give rise to a finite magnetisation, namely 
/ dxxp^{x) = 0. At the critical coupling Qc such solution will get unstable in favour of 
a non-symmetric one. In order to tackle this problem we first solve for the symmetric 
stationary solution employing a Fourier series expansion 

+00 

P*i^) = ^ 51 '^fe exp{inkx) . (7) 

V ^ k=-oo 

Here, because of normalisation cq = 1/ -\/2, and all the coefficients Ck = C-k are real 
because of symmetry. With the abbreviation 



Fkk'ig) = J exp[in{k'x - kgf{x))]dx 



4 cos (2vrF/3) ^^K^' + ?,kg)/?,] ^sin [vr(fc- - 3fc,)/3] 

^ ' ' n{k' + 3kg) 7i{k' - 3kg) ^ ^ 



the fixed point equation reads 




Ck = 7^ ( J2 ^kk'i^ - ^9)ck' I ( J2 Fkk"{.g)ck" . (9) 

k" J 
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Although such an equation can be numerically solved by iterative methods, it is 
convenient to look for an approximate analytical solution. As Fkk' + -Pfc_fc' = holds 
whenever the index k' is not an integer multiple of three, the right hand side of equation 
(y) contains only Fourier coefficients of the form Ck'=3i for symmetric densities. If we 
truncate the expansion at the fifth mode, then the right hand side of eq.(^ only 
contains the coefficient C3. The latter is self consistently determined by 

V2C3 = ^ (F3o(1 - 3g) + [Fssil - 3g) + ^3-3(1 - 3g)]V2c^) 

[fM + [FM + F3_3(^?)] V2C3)' . (10) 

The polynomial (10) admits two real solutions in the whole range < g < 1/3, but only 
one of them is smaller than 1/2. The remaining Fourier coefficients up to order five are 
now obtained if we plug in the expansion (^ with the solution of eq.(lO) into the right 
hand side of eq.(^. The analytical expression thus obtained coincides up to five digits 
with the full numerical solution of eq.(^ if many modes are taken into account. 

We are now going to evaluate explicitly the critical point. It corresponds to the 
value of coupling where the symmetric solution loses its stability. Considering therefore 
small deviations from it Pn = P* + Spn, and expanding eq.(|D we obtain 

6pn+l{x) = {C6pn){x) + {C[6pn,Spn]){x) + {V[6pn,6pn,6pn]){x)-\ .(11) 

The deviations obey the constraint / dx 6pn = because of the normalisation condition. 
The stability properties are determined by the eigenvalue problem for the linear operator 

{C6p){x) = J dyA{x,y)6p{y) . (12) 
Its kernel reads 

A{x, y) = J dzidz2dzs {5 [x - T{y, Zi, Z2, z^)] + 35 [x - T{zi, y, Z2, Z3)]) 

P*{Z1)P*{Z2)P*{Z3) . (13) 

Under quite mild conditions on p* (e. g. continuity is sufficient) the kernel is continuous. 
Hence the operator (p^ is compact and its spectrum consists of isolated eigenvalues 
which accumulate at most at zero 0. In addition, the bifurcation behaviour of the 
symmetric density is essentially identical to bifurcations in low-dimensional dynamical 
systems. In particular, the instability is typically caused by a single eigenvalue A crossing 
the unit circle in the complex plane. 

The numerical solution of the eigenvalue problem is accomplished by representing A 
on a truncated Fourier basis, and diagonalising the resulting finite-dimensional matrix. 
The spectrum thus obtained at the critical coupling is displayed in figure |I]. An isolated 
eigenvalue attains the unit circle along the real axis and the corresponding critical 
eigenfunction v{x) is odd with respect to the space inversion v{x) = —v{—x) (cf. figure 
^ f. Notice that, as the kernel (13) is not symmetric and the corresponding linear 

f Since the full equation (^) is homogeneous of order four the linear operator ( p^ ) admits a Goldstone- 
like mode with eigenvalue A = 4 and eigenfunction p^(x). 
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operator (|T2|) is in general not selfadjoint, the critical left-eigenfunction w{x) may differ 
from v{x). From the numerical diagonalisation of the matrix the critical coupling is 
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Figure 1. Spectrum of the operator ([l^) at the critical coupling gc- 
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Figure 2. The critical eigenfunction v{x) corresponding to A = 1. 
very accurately determined to be gc = 0.1496 . . ., which remarkably deviates only 6% 
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from the value 0.1584. . . of the critical point obtained in numerical simulations of the 
full map lattice with asynchronous updating [Q. This coincidence is consistent with the 
idea that the latter dynamic rule induces less spatial correlations than the synchronous 
one. 



3. Bifurcation analysis 

The bifurcation scenario described at the end of the preceding section is of the pitchfork 
type. Hence the growth of the critical mode beyond the instability is expected to scale as 
(g — gcY^"^. We stress that such a critical behaviour follows solely from the compactness 
of the linear operator. Therefore, we expect that every mean-field approximation of the 
Perron-Frobenius equation shares the same critical exponents. 

To support such general remarks we are now going to perform a complete normal 
form reduction of the full mean-field equation. For that purpose one needs the quadratic 
and cubic contributions of eq. ([llD which are given by 



{C[ip,'^]){x) = J dyidy2T{x,yi,y2)(p{yi)^{y2) (14) 

{V[(p,^p,uj]){x) = J dyidy2dy3A{x,yi,y2,y3)viyi)ijiy2)ujiy3) ■ (15) 
The corresponding kernels read for example 

r(a;,2/i,?/2) = J dzidz2[35{x -T{yi,y2,zi,Z2)) + 35{x -T{zi,yi,y2,Z2))] 

p^{zi)p.^{z2) (16) 
A{x,yi,y2,y3) = J dz[6{x - T{z,yi,y2,y3)) + 36{x - T{yi, z,y2,y3))] P*{z) ■ (17) 

The centre manifold tangential to v{x) is expanded as 

Spnix) = anv{x) + Q;^r(x) + als{x) H , (18) 

where the scalar an denotes the coordinate on the one dimensional manifold. To remove 
the ambiguity on the transversal vectors r{x) and s{x), we chose {w\r) = 0, {w\s) = 
with respect to the canonical bilinear form (ipl^p) = J dx il){x)ip{x) . The time evolution 
on the centre manifold obeys 

a„+i = Aan + Bal + Cal + ■■■ . (19) 

If we now plug in eqs. (PTSp and (|TU|) into eq.(|TT]) and compare the different powers in a„, 
then the unknown expansion coefficients in eq.(|l9D will be fixed. To the first order we 
just obtain the eigenvalue equation 

{Cv){x) = Av{x) (20) 

so that A = \. To the second order we obtain a linear inhomogeneous equation 
determining the transversal vector r 

A^r{x) - (£r) (x) = {C[v, v]){x) - Bv{x) . (21) 
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Since the operator on the left hand side becomes singular at the critical point g = g^, 
but the solution has to stay regular, the right hand side must obey a Fredholm condition 
at g = gc- In particular, the right hand side has to be orthogonal to w{x). Moreover, 
as the first term on the right hand side drops from this condition by symmetry, we are 
left with B = 0. Finally, at the third order we obtain 

Ah{x) - iCs){x) = {C[v,r]){x) + {C[r,v]){x) + {V[v,v,v]){x) - Cv{x) . (22) 

The Fredholm condition ensuring a regular solution at criticality determines the cubic 
coefficient as 

C{w\v) = {w\C[v,r]+C[r,v]+V[v,v,v]) . (23) 



Eq.(p3D is easily evaluated using eqs.(14), (15), and (|2lD . If the representation in terms 
of Fourier modes is employed, then all matrix elements can be expressed in terms of 
(8). Of course the modulus of C has no special meaning, since it depends on the 
normalisation of v{x). We just chose / dxv'^{x) = 1 and obtain C = —8.889 x 10~'^. 
Hence, the pitchfork bifurcation is supercritical as already mentioned above. The 
stationary distribution in the supercritical region is thus readily evaluated from the 
normal form (|T^ as p*(a;) + y^(A — !)/(— C) v{x), where p* is the (symmetric) stationary 
density at the critical point gc- Accordingly, the critical mode and the magnetisation 
grow like (A — 1)^/^ ~ ig~gcY^'^ beyond criticality. In addition, it is important to notice 
that the system is close to a super-subcritical transition. In fact, if one considers the 
definition (^31) for arbitrary coupling g, then a change in the sign of C occurs slightly 
below the transition point gc aX g = 0.1477.... For that reason, the correct scaling 
behaviour occurs only in a narrow region beyond g^ and the evaluation of the critical 
exponents from a numerical solution of eq. (^) is almost impossible to perform. The last 
point again emphasises the importance of our analytical results. 

Sofar we have dealt with the critical behaviour of the order parameter only. Let 
us discuss now the counterpart of the static susceptibility in equilibrium systems. In 
analogy with the latter, we may tentatively define it as the derivative of the order 
parameter with respect to a suitable "symmetry breaking field". This amounts to 
introduce some external parameter h spoiling the symmetry of the original single site 
map. Although not unique, a natural choice is for example to replace / in eq.(|l]) with 

r-3a;-2 for -1 < x < -1/3 

fh{x) = I 3x for -1/3 < a; < 1/3 . (24) 

[-3{l + h)x + 2-h forl/3<x<l 

Here the transition rate from [0, 1] to [—1, 0] is lowered upon increasing the parameter h, 
thus mimicking the effect of a static magnetic field in the Ising system. The mean-field 
susceptibility can then be defined as 



d{x) d 

X ■ = 



o f 

— J dxxp^{x) , (25) 



h=Q '^"^ " h=Q 
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where p^,(x) denotes now the fixed point solution of eq.(|^) with the single site map (p^). 
If we suppose that such a solution depends on /i in a differentiable way, we can take the 
formal derivative of eq.(BD obtaining 



dh 



h=0 



x)= - dyodyidy2dys5'[x - Th=oiyo,yi,y2,y3)] 



dTh{yo,yi,y2,y3) 



dh 



p*{yo)p*{yi)p*{y2)p*{y3) , (26) 



h=0 



where the definition of ([T2|) has been also taken into account. In order to solve this 
equation for dp^/dh\h-o: we expand both the latter quantity and the the right hand 
side of eq.(26) in terms of eigenmodes of C. As only one eigenvalue crosses the unit 
circle and the right hand side of eq.(26) remains bounded under quite mild conditions 
on p^,, all the coefficients of such expansion remain bounded in the vicinity of the critical 
point, except for the one of the critical mode v{x). This coefficient develops a singularity 
of the form |A — which carries over to dp^:/dh\h=o, provided that the expansion in 
terms of eigenmodes converges absolutely. 

The same result can of course be derived along the lines of bifurcation theory. Since 
the instability without the symmetry breaking field is governed by the pitchfork normal 
form and only a single eigenvalue becomes critical, one expects from the very beginning 
that the symmetry breaking unfolds the normal form to the cusp case (cf. 0). That 
indeed occurs if one follows the normal form reduction presented above. Altogether, we 
can conclude that the static susceptibility diverges as \g — gc\~^, in accordance with the 
simple mean-field theory for equilibrium systems. 



4. Conclusions 



We have shown that a mean-field approach for the Miller-Huse model reproduces 
the critical behaviour known from simple mean-field approximations in equilibrium 
statistical mechanics. The value of the critical exponents basically originates from the 
compactness property of the linearised operator. Since such a property holds quite 
generally, our result is supposed to apply for a large class of coupled map lattices and 
almost all kinds of mean-field approximations. Corrections to mean-field scaling can 
only come from a continuous spectrum of the full Perron-Frobenius operator, similar 
to findings in equilibrium statistical mechanics, e. g. for the two dimensional Ising 
model. The coincidence with equilibrium mean-field theories is far from being obvious. 
For instance, globally coupled maps behave quite different compared to the mean- 
field approximation of the type employed here (cf. and references therein). That 
observation is in striking contrast to equilibrium statistical mechanics, where mean- 
field approaches and long-range coupled models are often equivalent. 

As mentioned in the beginning, the critical behaviour of the full map lattice depends 
on the updating rule. In particular, the sole occurrence of equilibrium Ising exponents 
for asynchronous updating was attributed to a kind of destruction of coherence by the 
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updating rule. Of course, the plain mean-field approach cannot explain such differences 
in critical exponents at all. Nevertheless, the mentioned interpretation is fully consistent 
with the fact that the mean-field approach (which neglects all correlations) yields a good 
estimate of for the model with asynchronous updating. 

At equilibrium, the static susceptibility can be expressed in terms of the spatial 
correlation function by virtue of the properties of the canonical distribution. This is of 
course no longer true for out of equilibrium systems as the coupled map lattices. Hence, 
it would be tempting to check whether even in this case the critical behaviour of the 
susceptibility coincides with that of spatial correlations. Such a coincidence, which of 
course requires quite accurate numerics, would indicate a relation between response and 
correlations in nonequilibrium systems too. 
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